Current-induced torques in continuous antiferromagnetic textures 
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We study the influence of an electric current on a continuous non-collinear antiferromagnetic 
texture. Despite the lack of a net magnetic moment we find that the exchange interaction between 
conduction electrons and local magnetization generally results in current-induced torques that are 
similar in phenomenology to spin transfer torques in ferromagnets. We present the generalization 
of the non-linear sigma model equation of motion for the Neel vector that includes these current- 
induced torques, and briefly discuss the resulting current-induced antiferromagnetic domain wall 
motion and spin-wave Doppler shift. We give an interpretation of our results using a unifying 
picture of current-induced torques in ferromagnets and antiferromagnets in which they are viewed 
as due to the current-induced spin polarization resulting from an effective spin-orbit coupling. 
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Introduction — Magnctorcsistive phenomena in conduct- 
ing ferromagnets, resulting from the interplay between 
spins of conduction electrons and the magnetization, are 
well known. They play a major role in new technologies 
and are key to the recently awarded Nobel prize for Gi- 
ant Magnetoresistance[lJ,l2|. Recently, the effect of a spin 
current on magnetization dynamics — called spin trans- 
fer — has been investigated in detail 0,0, HLH ■ Although 
applications, mainly to memory storage technology, are 
an important driving factor behind this research, spin 
transfer is also physically interesting in its own right. It 
can be understood as follows. 

Consider a conducting ferromagnet far below its criti- 
cal temperature such that it is described by a unit vector 
fl(x, t) in the direction of magnetization. Its dynamics in 
the presence of an effective field, H c g, is determined by 



by contributing to the effective field as 



— = n x h cS - ,m x v 2 n, 

at 



(1) 



where we have ignored magnetization relaxation and the 
effective field contains in first instance contributions from 
the anisotropy of the system and the external field. Fur- 
thermore, J is the exchange constant favoring alignment 
of neighboring spins. In the presence of conduction elec- 
trons the effective field contains an additional contribu- 
tion due to the s — d exchange coupling of the magnetiza- 
tion with the spin density s of the conduction electrons 
given by — J dxAs- O/a 3 , with A the exchange splitting 
and a the lattice constant. (Although the s — d model 
is convenient to illustrate the physics involved the con- 
clusions drawn are qualitatively valid for other models as 
well.) In equilibrium, i.e., without a current, this spin 
density will be in the plane of the magnetization. How- 
ever, when a current is applied, the spin density acquires 
a component (s±) perpendicular to the magnetization di- 
rection f2. It is this component of the spin density that 
leads to current-driven dynamics of the magnetization Q 



Aa 3 



(2) 



To lowest order in the current and the gradient of the 
magnetization, and using an adiabatic assumption based 
on the fact that the electron dynamics is fast compared to 
the time variation of the magnetization, the out-of-planc 
spin density is 



Aa 3 



fix [v s ■ V fi 



(3) 



where the velocity v s , defined via the above equation, is 
proportional to the electric current. Inserting the above 
contribution to H c g in Eq. ([T]) wc find that in the steady- 
state transport situation the contribution of the current- 
induced torques to the equation of motion for the mag- 
netization direction is given by 



dSl(x, t) 



at 



V)fl(x,t). 



(4) 



The current-dependent velocity is in the absence of spin- 
orbit coupling and spin-flip scattering fixed by spin 
conservation, as follows. Integrating Eq.(|3]) over the 
length of the ferromagnet in the direction f of the cur- 
rent leads to the total change in angular momentum 



Jdrdn(x,t)/dt 



current 

current flows from in to out 



fi 



- fi 



where the 



By spin conservation, 
this change in angular momentum is proportional to 
the change in the spin current Jf after it has passed 

through the ferromagnet, i.e., fdrd£l a (x,t)/dt = 

current 

J£ in - J£ out . Using that Jf = PJ c fi a , with P the po- 
larization of the charge current J c in the ferromagnet, 
we find that v s oc PJ C . This argument shows that the 
current-induced torque in Eq. (j4|) results from transfer of 
angular momentum from conduction electrons to magne- 
tization, and is hence called a spin transfer torque. 
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Well-known examples of dynamics resulting from these 
spin transfer torques involve spin waves and domain 

waiis d i, EI EI Ei El EI El EI Ei El, El- For 

spin waves a Doppler shift in the dispersion relation is 
found once a current is applied [1, Q . This shift is lin- 
early proportional to k where the proportionality con- 
stant is given by v s and the dispersion that follows from 
Eqs. ([T][3])is then given by hio = Jk 2 +v s ■ k, where we 
ignored anisotropy and external fields. Recently exper- 
iments have successfully measured such current-induced 
spin-wave Doppler shifts 19 . Current-induced domain 

3, Ei, EI El, EI Ei is un- 
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wall motion [ 

derstood as follows. In the absence of any pinning and 
damping, Q,o(x — v s t) is a solution to Eq.Q, where fto(x) 
is a timc-independent solution of Eq.|T|). For the case 
that flo(x) corresponds to a domain wall texture this 
means that the velocity of the domain wall is just v s , 
an estimate which turns out to be reasonable even when 
pinning and damping are present H 11 1. These results 
are also understood by realizing that Eq.(j4]) shows that 
it is possible to get the equation of motion of O with an 
applied current by substituting — ► ^ + v s ■ V into the 
equation of motion without current. 

In antiferromagnets adjacent lattices sites have op- 
posing magnetic moments and thus form two sublat- 
tices with opposite magnetization leaving no net mag- 
netization. They are characterized by a Neel vector 
fij = ( — lytlj where j labels the lattice sites. Its equa- 
tion of motion is a non-linear sigma model in the presence 
of an effective field, i? c ff , and is given by 



on 



d_ 

di 



(n x H cff j 



c 2 V 2 n 



0, 



(5) 



to lowest order in H e e and where c is the spin- wave ve- 
locity. In this Letter we show that in the presence of an 
electric current in an antiferromagnetic conductor there 
is, to lowest order in Neel-vector gradient and the cur- 
rent, an out-of-planc spin density given by 



<*0 = 



Aa 3 



T n x 



(v ■ V) fi, 



(6) 



where v, defined with the above equation, is proportional 
to the current and the antiferromagnetic equivalent of 
the velocity v s that was introduced in describing current- 
induced torques in ferromagnets. Although the form of 
the spin density is the same as for a ferromagnet, the 
velocity v cannot be determined from macroscopic spin 
conservation arguments, as in the case of the ferromag- 
net, but instead needs to be determined by a microscopic 
calculation (such a calculation is discussed in detail be- 
low). Using the equation of motion, Eq.([5]), and the spin 
density in Eq. © as a contribution to _ff c ff gives the equa- 
tion of motion 



d 2 n 







+ ^-(v- V)n- c 2 V 2 r 



dV 







(7) 



where we ignored anisotropy. The antiferromagnetic spin 
wave dispersion resulting from this equation is, again 
ignoring anisotropy and external fields, given by ftio — 
ck + v ■ k/2 to linear order in current. This shows that 
the dispersion has a similar Doppler shift as the ferromag- 
netic case 8, 9|. Similarly like in the ferromagnetic case, 
the equation of motion in Eq. ([7]) allows for co-moving 
solutions no(x — vt/2) (to first order in v), from which we 
conclude that antiferromagnetic domain walls [Iflj move 
with velocity v/2 proportional to the current in the ab- 
sence of pinning and damping. Note that the above equa- 



tion of motion is obtained by replacing 



dt 



■V/2 



in Eq. ([5]) and keeping terms to first order in v. 

So far we have suggested a form for the spin density of 
the conduction electrons in an antiferromagnetic metal in 

the presence of current, i.e., (s_j_) cx ft x yv ■ V^j n. In the 
remainder of this Letter we prove this result and show 
that v is generally not zero and can be of the same order 
of magnitude as v" s , which is of the order of 1 — 100 m/s for 
a typical ferromagnetic alloy like permalloy. We calculate 
v for a toy-model antiferromagnetic metal. Further, we 
explain our results using a unifying picture of current- 
induced torques in ferromagnetic and antiferromagnetic 
textures, that is built on the notion of effective spin-orbit 
coupling induced by a non-collinear magnetic texture. 

Previous work on current-induced torques in antifer- 
romagnetic metals considered mainly single-domain lay- 
ered structures [U, 22, 2l 24, (2{| and the situation of 
an antiferromagnetic domain-wall was considered from 
an ab initio point-of-view [26j |. In this Letter we derive 
the general phenomenology of current-induced torques 
in spatially smooth and slowly- varying antiferromagnetic 
textures. 

Toy-Model Antiferromagnetic Metal — To confirm the 
conjecture (s±) cx nx (v- V)n and determine the current- 
dependent velocity v we consider a toy model of a con- 
ducting antiferromagnet within the Green's function for- 
mulation of Landauer-Buttiker transport theory 27J • The 
tight-binding hamiltonian used in this theory is given by 
H = H s + H L + Hi where 



He 



E 

<j,j'>;<r 



(8) 

is the system Hamiltonian with t the nearest-neighbor 
hopping amplitude and A the exchange energy, and ^ >cr 
are tpj a are the electron annihilation and creation op- 
erators. Hl and Hi are the Hamiltonians for the leads 
and for the coupling between the leads and the system 
respectively and are similar but with A = 0. To realize a 
transport current I these leads have a chemical potential 
difference of e\V\. The magnetization texture is set to 
tij = ((-l) J sin(2^),0,(-l)Jcos(^)) where A is the 
wavelength of the antiferromagnetic texture. See Fig. [1] 
for an illustration. This magnetization texture is in the 
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FIG. 1: Illustration of a smooth antiferromagnetic magneti- 
zation texture. 



x — z plane and in equilibrium, without current, the sys- 
tem only has nonzero spin densities s — (i/;J.T (Tj(J ''f/v) in 
this plane. For nonzero voltage we find a spin density s± 
in the y-direction that is constant in position, in agree- 
ment with Eq. (J6j> . According to Eq. ([6]) the velocity, 



v/n 



AAJgx) 
2irh ' 



(9) 



should be independent of A in the long- wavelength limit. 
Fig.[2]confirms this. The fluctuations for small A in Fig. [2] 
are due to finite size effects. The ferromagnetic result is 
also shown in this figure. Note that the long-wavelength 
limit is reached for smaller A in the ferromagnetic case, 
compared to the antiferromagnetic one. This is under- 
stood as the unit cell of the antiferromagnetic system is 
effectively doubled because of the opposing magnetiza- 
tion on neighboring sites. 



AASi /la 
2irfi 1 \e\ 




FIG. 2: Plot of the spin density versus the texture wavelength, 
A. The solid line corresponds to the antiferromagnet and the 
dashed line represents the ferromagnet. 



Now that we have established that v approaches a 
constant in the long- wavelength limit, we take its long- 
wavelength limiting value as its definition and study its 
dependence on A. The result is shown in Fig. [3] ft is 
clear from this figure that v depends strongly on the value 
of the exchange constant. This dependence reflects the 
strong dependence of the quasi-particle band structure 
on A. The ferromagnetic case is also shown for compar- 
ison. Note that in the ferromagnetic case the velocity 
for increasing A quickly obtains its value determined by 
spin conservation. In the antiferromagnetic case the ve- 
locity is not determined by a macroscopic conservation 




FIG. 3: Plot of the velocity, v, versus the ratio of the exchange 
constant and Fermi energy, — . The solid line represents the 
antiferromagnet and the dashed line represents the ferromag- 
net. 



law. Note, however, that the velocities v s and v, respec- 
tively parametrizing current-induced torques in ferro and 
antifcrromagnets, are of the same order of magnitude. 
Effective Spin- Orbit Coupling — The above numerical 
analysis confirms that when a current is applied to an 
antiferromagnetic metal there is a nonzero out-of-planc 
spin density that drives current-induced magnetization 
dynamics. We now give an interpretation for this spin 
density in terms of a current-induced spin polarization 
resulting from an effective spin-orbit coupling. We start 
with a system of electrons (mass m), moving in an anti- 
ferromagnetic texture n(xj) and scalar potential V(£j), 
described by the Hamiltonian 

H=^- + V(x' J )-^(-iyn(x' J )-T, (10) 

with Xj = ja the position of the j-th lattice site. We align 
the spin quantization axis of the conduction electrons to 
the local Neel vector by applying a SU (2) transformation 
— ► R\ip) to the wavefunction with ■ fR = t z 

that therefore diagonalizes the spin part of the Hamilto- 
nian and gives the effective Hamiltonian 



H eS = H Q -2iJ»a»(Z), 



(11) 



to first order in the gradient of ft, where Hq is the Hamil- 



hdkl with e fe 



tonian in Eq. (|10[) with n — z. Here, Jjf a = ^t^Vq 

the spin current with the velocity uo, Q 
the dispersion corresponding to H and the gauge fields 
a£(x) = in x V^n| r [28j|. Substituting the latter into 
Eq.(8) leads to the effective Hamiltonian 



H 
H 









x 


dk a a 






J 



B e s ■ T, 



112) 



where B c tf is a momentum-dependent fictitious magnetic 
field that can be viewed as an effective spin-orbit cou- 
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pling resulting from the non-collinear antifcrromagnctic 
texture. (Note that this effective spin-orbit coupling is 
different from the result of Ref. [29[ for a collinear antifer- 
romagnet.) Since the effective magnetic field is linearly 
proportional to the velocity operator, when an electric 
field is applied B c fi is nonzero. The resulting Zeeman 
splitting of the electron spins due to B e g results in a 
current-induced spin polarization that is aligned with the 
effective magnetic field and, since i? c ff oc n x Vn conse- 
quently is perpendicular to both the magnetization and 
its gradient, ft is this current-induced spin polarization 
that contributes to the effective field for the magnetiza- 
tion and leads to current-induced torques. The above ar- 
gument holds in an analogous form for ferromagnets and 
confirms that the magnetization-direction dependence of 
(sj_) for ferromagnets, is the same as its Neel vector de- 
pendence for antiferromagnets. Note that at the level of 
the effective hamiltonian in Eq. (p~2|) the main difference 
between antiferromagnetic and ferromagnetic case is the 
difference in Hq, which in the ferromagnetic case con- 
tains a constant exchange splitting and in the antiferro- 
magnetic case an alternating one. Current-induced spin 
polarization has been studied in paramagnetic semicon- 
ductors l3ll I32I ] , which shows that they do not require 
a net nonzero exchange splitting, which, in turn, explains 
why the current-induced torques in antiferromagnets are 
generally nonzero. 

Discussion & Conclusions — So far we have ignored 
conduction-electron spin relaxation which in the ferro- 
magnetic case is known to result in an additional contri- 
bution — /3VO to the spin density that is parameterized 
by the dimensionless constant (3 [33L 34 1. As we have 



shown the spin density has, without spin relaxation, the 
same form for both ferromagnetic and antifcrromagnctic 
textures. It is therefore reasonable to assume that spin 
relaxation will lead to a similar correction for antiferro- 
magnets. Hence, we expect that the spin density takes 
the form (sj_) oc fi x (jj ■ V^j n + /3 a fm ■ V^j n when we 
include conduction-electron spin relaxation that leads to 
the correction phcnomcnologically parameterized by the 
dimensionless constant /3 a f m . This will change the equa- 
tion of motion Eq.((7|) which now becomes 



d 2 ri d 
~dt 2 ~ + d~t 



d f 

c 2 V 2 n+/3 afm — {n x [v ■ V]n 



= 0. 
(13) 

In future work we intend to explore the consequences of 
this equation of motion for the current-driven motion of 
antifcrromagnctic domain walls, taking into account also 
anisotropy and magnetization relaxation. Further studies 
will also include exploring other aspects of the effective 
spin-orbit coupling discussed above. 

In conclusion, we have derived a general equation of 
motion for the influence of an electric current on the 
long- wavelength dynamics of a smooth antifcrromagnctic 
magnetization texture. Although we have only presented 



numerical calculations for a simple model, we believe that 
the form of this equation of motion is quite general, be- 
cause the arguments leading to Eqs. (fT^|) and fT3")) do not 
depend on details. 
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